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Abstract 



Q I We study static solutions of the Skyrme model on the two-sphere of radius L, for 

r"! I various choices of potential. The high-density Skyrmion phase corresponds to the 

' , ratio j3 = L/(size of Skyrmion) being small, whereas the low-density phase corre- 

D , sponds to p being large. The transition between these two phases, and in particular 



the behaviour of a relevant order parameter, is examined. 



X, 

1 Introduction 



Topological solitons are usually studied on flat space, but there are several reasons 
why the curved-space setting is interesting (apart from the most obvious one of pos- 
sible cosmological applications). First, it often reveals useful mathematical features; 
there are now two length scales, namely the size of the solitons and the radius of 
curvature of the underlying geometry, and there is an interplay between these scales. 
For monopole-like objects such as as vortices and Yang-Mills-Higgs monopoles, one 
can work on hyperbolic space: in the case of vortices, the equations for static vor- 
tices become integrable for a particular value of the curvature [||], |||]; and for YMH 
monopoles, the system is simpler on hyperbolic space than on flat space [^. For 
textures, such as the Skyrmions discussed in this paper, one may study the problem 
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on either hyperbolic spaces or spheres (ie on spaces of constant negative or positive 
curvature), and this similarly throws up interesting geometrical features Q, j^, @, 

0, i, i, 0. 

There are also more physical reasons for studying textures on compact spaces, 
namely that this models a finite density of solitons, and the transistion between 
the high-density and the low-density phases. One may do this without introducing 
curvature, by taking space to be a flat torus; but it is in some ways simpler to use 
a sphere. This was first done for the Skyrme model in three spatial dimensions, 
which is a model of nuclear matter. But there has also been interest in the two- 
dimensional version, both as a 'toy' model and also because it has several potential 
applications in condensed-matter physics. In this paper, we study the Skyrme model 
on the two-sphere of radius L, and in particular the dependence on the ratio /? = 
-L/(size of Skyrmion). 

In the three-dimensional Skyrme model (without a pion mass term) , the Skyrmion 
for small (3 is homogeneously spread out over space; when f3 reaches a certain critical 
value, the Skyrmion localizes around a point ||5[ , |]lO| . If the winding number is greater 
than one, or if a zeroth-order potential term (such as a mass term) is present, then 
this transistion is less clear-cut. In the two-dimensional case, such a potential term 
is obligatory, and so this smearing effect is inevitable (c/ Q). The quantity $ which 
monitors the deviation from homogeneity tends to zero as f3 does, but ^ is never quite 
zero, and so is not an order parameter in the strict sense. There is, however, another 
quantity ^ which does serve as an order parameter, namely one which monitors 
the breakdown of a certain reflection symmetry [10|; and the associated phase 
transition is sharp (as long as the potential has this symmetry, which was not the 
case for the system studied in [^). We shall investigate these phenomena, for several 
different potentials, by finding rotationally-symmetric solutions numerically. 



2 Skyrme models on the two-sphere 

In this section, we review the Skyrme system on the two-sphere, the topological 
Bogomol'nyi bound on its energy, and the imposition of rotational symmetry. Let M 
denote the standard 2-sphere of radius L, with local coordinates {x^^x^) and metric 
ds^ = Qjk dx^ dx^. The area element on M is w = h dx^ Adx"^ , where h = [det{gjk)]^^'^ ■ 
The Skyrme field (j)is a map from M to the unit sphere S'^, and may be thought of as a 
unit vector field (j){x^) (so 0-</> = 1). We are considering only the static problem, so (p 
depends only on the spatial variables x-', and not on time. Let dj(j) denote the partial 
derivative d(p/dx^. The winding number (degree) of cp is denoted by the integer k. 
If we define the function p[(j)] as the triple scalar product p[(j)] = 2(j) ■ (9i</>) x (920), 
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then (assuming that the cover almost all of M) we have 



k = ^ [ p[(f]dx^ dx^. (1) 



Sit 



The (normalized) energy Ef^ of (p is taken to be a sum of three terms Ej^ = E^'^'^ + 
£;{4) -|-^(o) ^ with each being the integral over 5^ of the corresponding density function: 

for n = 2,4,0. The densities £^^\(p) are given by 

= ^g'HdjMdk^), (3) 

£:(^)[<^1 = lih-'p', (4) 

= aVi$), (5) 

where 7 and a are positive constants, and where V{(j)) is some suitable function of 
(f) (not involving its derivatives). We regard the normalized energy as dimensionless. 
The constants 7 and a have dimensions of length^ and length"^ respectively; the 
dimensionless combination 07 determines the energy of a soliton, while its size is 
of order (7/a)^/^. In flat space, this latter quantity simply sets the length scale; 
but on the sphere M, there already is a length scale L. So one has a dimensionless 
parameter /3 = L(a/7)^/^, which is the ratio between the size of space and the size 
of the soliton. From now on, we shall set 7 = 1. 

The energy of any configuration of nonzero degree satisfies a generalized Bo- 
gomol'nyi bound |11|. First, E^'^^ satisfies the standard 0(3) sigma-model bound 
£^(2) > 1. To get a bound on E^^'^ + E^^\ write f^^) = and observe that 

> — — T / pG dx^ dx"^ 

= ^ [ G{$) ^-d^xdf, 
4vr Js2 

the last equality follows from the degree theorem J 4>*{G0,) = k J GO, with = 
(j) ■ d(f) X d(j) being the area element on the target sphere S^. If G depends only on 
the third component (p^ of 0, which is the case for all the examples that we consider, 
then the bound on the total energy takes the form 



1 

= 1 + j ^Vvi^ 



(6) 
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In this paper, we study rotationally-symmetric configurations. Let us use polar 
coordinates {0,{p) on M, so that ds"^ = L?'{d9'^ + sin^ 6*^99^) and h = L^sin^. We 
impose rotational symmetry by taking (j) to have the form 

4>= (sin(/)cos(A;y?),sin(/)sin(/cv3),cos(/)), (7) 

where k is an integer and / = f{6). Assuming that this profile function f{9) satisfies 
the boundary conditions /(O) = tt and /(vr) = 0, the configuration (0) has winding 
number k. Its normalized energy is 



ifr + k'('^)\(^'^y + 2L'£^'\cosf) sin6d9, (8) 
\smt^ / \ L smo J 



where /' = df /dO and where we are (as mentioned above) taking £^^^ to be a function 
of (p"^ = cos / only. 

We are interested in the stationary points of the functional (|8|), and we find 
these numerically (with an accuracy better than 0.1%). Two independent methods 
were used: first, a direct minimization of (^) by a conjugate-gradient method; and 
secondly, solving the corresponding Euler-Lagrange equations. 

One particular 'trial' profile which is of interest is the one corresponding to a 
rotationally-symmetric conformal map (c/ f^, @, ||lO|). Such a profile can be used 
for any value of k, but we shall restrict here to the case A; = 1. The profile function 
f{9) is then given by 

cot(//2) = Atan(0/2), (9) 

where A is a positive constant. Notice that if A = 1, then f{9) = it — 9, and this 
corresponds to the identity map from M = to S^. This identity map is only 
a solution of the equations of motion if S^^'^ is constant (ie not depending on 0). 
Another way of putting this is that if (and only if) £^^^ is constant, then Skyrmions 
spread out to fill the whole of space homogeneously, in the sense that the energy 
density £ is constant. The configuration @ has £'*^^^(A) = 1 and £'*^^^(A) = (A^ + 
A2 + 1)/(6L2a2). Notice that E^'^^X) has its minimum at A = I; but the location 
(and existence) of a minimum in the total energy E{X) depends on the choice of the 
potential term E^^\ The question is: for a given L and choice of E^^\ what are the 
minima of E{X), and in particular are they at A = I or at A 7^ I? 

We now wish to define the two quantities <I> and ^ mentioned in the introduction. 
The first of these is 
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This is zero for the homogeneous configuration f{6) = tt—9, and so gives an indication 
of the deviation from homogeneity. The second quantity is 

2 

(11) 



a3\2 



1 



47rL2 



M 



This is zero for configurations which have the reflection symmetry (p^ 
nonzero if the Skyrmion is locahzed around the point = or 9 = n. 



but 



3 Two asymmetric examples 

In this section, we give a brief discussion of two systems, arising from two possible 
choices of S^'^h Each is 'asymmetric', in the sense that neither posesses the symmetry 
(j)^ ^ —(j)^. So a Skyrmion in these systems wih never have ^' = 0; the field prefers to 
be near = (0, 0, 1) rather than (p = (0, 0,-1). The first example (for corresponding 



flat-space studies see [|^, [|T3[) is motivated by the question: can one saturate the 
Bogomol'nyi bound (^)? We shall not give a complete analysis of this question here; 
but in the k = 1 case, it is easy to derive the following fact: the configuration (^) 
saturates the Bogomol'nyi bound if and only if is given by 

■(A2 + i)_(a2_i)0" 



r(o) _ i 
^ " 2 



2LX 



(12) 



(Notice that the A = 1 case is degenerate: £^^^ is then constant.) For this system, 
the bound is 

^> 1 + (A^ + A2 + 1)/(3L2a2). (13) 
If, for example, we choose A = 2L, then the configuration (|^ is a solution with energy 



in the limit L ^ oo, this corresponds to fl^, [|l^] with their parameters 9i and 92 
set equal to 0.5 and 1 respectively (note also that their expression for energy E is 27r 
times ours). The energy ( p^ is a monotonic-decreasing function of L; this remains 
true for variations of (|l2|) such as S^'^^ = (1 — (f>^)^, and for k > 1. 

In this system, there is a repulsive force between solitons |12|, [|^], |14]. Con- 



sequently, one expects that the minimum-energy configurations for k > 1 will not 
be rotationally-symmetric (or, equivalently, that all rotationally-symmetric k > 1 
configurations are unstable). So investigation of the k > 1 sectors requires a full 
two-dimensional study. 
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The comments of the previous paragraph apply equally well to our other example 
in this section, which was the one studied in |jl^, [|^] (on flat space) and P 
(on the sphere). For this system, we have V{cl)^) = 1 — 0^. The Bogomol'nyi bound 
(|^ becomes E/. > 1 + 4-^/a/3 (note that this is stronger than the bound given in 
Q). Here we give a brief description of the k = 1 case, taking a = 0.1. For the 
A-approximation iM), we have 



^W(A) 



2L^a 



1 + 



A2log(A' 
1- A2 



2^ 



(15) 



which is a decreasing function of A. The total energy -E'(A) has a minimum at A = 
Amin> with Amin > 1 depending on the parameter /? = La^^^. In figure 1, we plot 
the energy E{Xiam), the actual Skyrmion energy E = Ei (the minimum of (§)), 
and the quantities ^, <I> and Amiin all as functions of L. Notice that the lowest 
value of E (when L = 1.65 and /3 = 0.93) is 1.424, which less than 0.2% greater 
than the Bogomol'nyi bound 1.4216. At this value of L, we have Amin = 1-39 and 
-£'(Amin) = 1.425. As L ^ cxD, the energy of the Skyrmion tends to its flat-space 
value of 1.549 0. 

As was noted in Q, there is a regime (/? < 1) where the Skyrmion is approximately 
spread out over M and well-approximated by the configuration and one (/3 > 1) 
where the Skyrmion is localized and not well-approximated by (^. But the transition 
between these two regimes is not sharply-defined, and is not a phase transistion in 
the usual sense. This is partly due to the asymmetry of the system. In our next 
examples, we shall see a sharper transition. 



4 Two symmetric examples 

In this section, we study the systems which are defined, respectively, by V{(j)^) = 
1 — (^^)^ and V{4>^) = (</>^)2[l — (</>^)^]. In these systems, Skyrmions attract one 
another |jl8|, |lj]; and consequently the solutions with k > 1 have rotational sym- 
metry. Since these potentials have the symmetry cj)^ i— > —0^, one expects (in the 
rotationally-symmetric case, and for large L) that there will be a stable Skyrmion 
solution localized at a point (either 9 = or 6 = n); and also that there exists a 
solution which is symmetric under 9 ^ tt — 6 (and which may or may not be stable). 
For small L, however, only the symmetric solution exists. 

First take V{(f)'^) = 1 — (0^)^, with the value a = 0.08873 (this corresponds to the 



value used in the flat-space study of |18|, and allows a quantitative comparison with 



its results). The Bogomol'nyi bound (^ is > 1 + tt^/oJ^ = 1.3309. Let us look 
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first at the k = 1 sector. For the A- approximation (P), we have 



ii;(°)(A) 



(1-A2)2 



2 



1 + A 
1-A2 



log(A^) 



(16) 



which is a positive function with a maximum at A = 1 (recall that £^(^)(A) has a 
minimum at A = 1). Notice that E{X~^) = E{\); this corresponds to the (j)^ i-^ — 0^ 
symmetry of E[(p] in this case. So A = 1 gives either a minimum or a local maximum 
of E{X), depending on the value of L. A simple calculation reveals that if < 5/2, 
then 

'^("-'-^-^ (") 

is a minimum of E{X); while if > 5/2, then there is a local maximum at A = 1, 
and the minima are at Amin > 1 and 1/Amin- Consequently, if we restrict to the 
special profiles (^), then there is a sharp phase transition at = 5/2 between a 
homogeneous phase (A = 1) and one where the Skyrmion localizes at a point (in this 
case, one of the poles) of the sphere M. Notice also that the lowest value of -E'(l), 
namely ^ = 1 + = 1.344, is attained when /J^ = 3/4 44> L = 1.705. 

This is just an approximation; the true Skyrmion energy is strictly less than 
E{X). However, there is indeed a minimum at L = 1.705 (with E = 1.341, only 
0.2% less than the A-approximation); and there is a transition at P'^ = 5/2 <^ L = 
2.304, as one can see from figure 2. Figure 2(a) depicts the energies of two static 
k = 1 solutions, as functions of L. The solid curve is the energy of the 'symmetric' 
solution, ie the one satisfying /{n — 6) = n — f{6). This is obtained by starting 
with the configuration f{9) = tt — 9 and iterating towards a stationary point; our 
procedures are able to converge to a saddle-point (which this is for large L). The 
dashed curve is the energy of an asymmetric solution, obtained by starting with 
a highly- asymmetric configuration. We see that for < 5/2, the two solutions 
coincide; but for > 5/2, the former becomes unstable (numerical experiment 
indicates that it is a saddle-point rather than a local minimum). The corresponding 
parameters ^ and $ for the minimum-energy solution are plotted in figures 2(b) and 
2(c); note that when < 5/2, we have = (the solution is symmetric), but $ 7^ 
(it is not homogeneous). There is a second-order phase transition at P'^ = 5/2. In 
figure 2(d), the Skyrmion energies (minima of E) are plotted for k = 1,2, 3. It might 
be noted that E^ gets very close to the Bogomol'nyi bound, at L = 1.7, L = 2.9 and 
L = 4.3 for k = 1, 2, 3 respectively (only 0.1% above for k = 2). As L ^ 00, we have 

El 

El 1.564, E2 1.405, ^3 ^ 1.371. (18) 
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Now let us turn to the system with V{((>^) = — {<t>^)'^], taking a = 0.04 

(which is consistent with the choice made in The Bogomol'nyi bound is E > 

1 + -v/2a/3 = 1.0943. In the k = 1 sector, the A- approximation (P) has 

AL'^aX^ [(10 + 18A2 - ISA^ - lOA^) + (3 + 21A2 + 2l\'^ + ^X^) \og{X^)\ 

^ = 3(A2 - 1)5 ' 

(19) 

which has a local minimum at A = 1, with S(°)(l) = 2L'^a/lb. Consequently, ^(A) 
always has a local minimum at A = 1. Notice that for L = [15/(4a)]i/'^ = 3.11 
(where the true energy has its lowest value 1.1026), the approximation has energy 
E{1) = 1 + Y^4a/15 = 1.1033. 

There are two values Li and L2 which are critical in the following sense. For 
L < Li, E{X) has only one stationary point (the minimum at A = 1). For L > Li, 
E(X) has three local minima: a symmetric one (A = 1) and a degenerate asymmetric 
one (A 7^ 1). For Li < L < L2, the asymmetric minimum has higher energy than the 
symmetric one: see figure 3(a), where E is plotted as a function of = (A^ - 1) /(A^ + 
1), for L = 7.4. Finally, for L > L2, the asymmetric minimum has the lower energy 
(c/. figure 3(b), for L = 8.4); and the energy of the symmetric minimum tends to 
infinity as L ^ cxd. 

This approximate picture suggests that for L large enough, there are two stable 
Skyrmion solutions (symmetric and asymmetric); but that only the latter survives 
in the limit as L ^ 00. Our numerical investigation shows that this is indeed the 
case. Results from minimizing (^ are presented in figure 4. In figure 4(a), the 
solid curve is the energy of the symmetric k = 1 solution; while the dashed curve 
is the energy of the asymmetric solution, which exists only for L > Li ~ 6. For 
Li < L < L2 ^ 6.9, the energy of the asymmetric solution is greater than that of 
the symmetric solution. Numerical evidence indicates that each of these solutions is 
stable (fe is a local minimum of the energy), at least in the rotationally-symmetric 
class. The stability of the symmetric solution, even for large L, is a consequence of 
the form of the potential V = (0^)^[1 — (<;/>^)^], which favours </>^ = 44> f{6) = 7r/2. 
It also favours (j)'^ = ±1, which stablizes the asymmetric solution. Between these 
local minima there will be a saddle-point solution, but we have not investigated this. 

Figures 4(b) and 4(c) show the quantities ^ and <I> for the minimal-energy k = 1 
solution; and in figure 4(d), the minimal energy is plotted for /c = 1, 2, 3. As L ^ 00, 
we have |14] 



El 1.281, E2 1.166, E3 1.143. (20) 
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5 Concluding Remarks 



We have studied rotationally-symmetric static solutions of the Skyrme model on the 
two-sphere of radius L. Even in the sector with winding number k equal to unity, the 
properties of the solutions, and of the transition between the high-density (small L) 
and low-density (large L) phases, depends crucially on the choice of potential term. 
It should be instructive to investigate the k > \ solution spaces in more detail 



(c/ fl^]). For symmetric potentials (such as those of section 4), this can 

be done within the rotationally-symmetric class; this is the analogue of the rational- 
map ansatz [IC] in three dimensions. But for asymmetric potentials such as those of 
section 3, one expects a different picture: the Skyrmions should separate and form 
some pattern on 5^ as their optimal configuration. Flat-space studies |jl^ indicate 
that this pattern will, in general, be rather complicated. 
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L L 

(c) $ (d) A„,i„ 

Figure 1: Properties of a single Skyrmion, for the system with V = 1 — (f)^, as a function 
of L. In (a), the dot-dashed Une is the Bogomol'nyi bound, the sohd curve is the actual 
energy, and the dashed curve is the best A-approximation. 
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(c) $ (d) Energies for fc = 1, 2, 3 

Figure 2: Skyrmion properties for the system with y = l — (^3)^, asa function of L. 
In (a), the dot-dashed hne is the Bogomol'nyi bound; the sohd curve is the energy of the 
symmetric solution; and the dashed curve is the energy of the asymmetric solution. Figures 
(b) and (c) plot the quantities \& and $ for the k = 1 minimum. Figure (d) gives minimum 
energies for /c = 1 (solid curve), k = 2 (dashed curve) and k — 3 (dot-dashed curve); the 
solid line is the Bogomol'nyi bound. 
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(a) Energy E{fi) for L = 7.4 (b) Energy E{fi) for L = 8.4 

Figure 3: Energy E{fi) of the approximation (|), for the system with V = (03)^[1 — {(psY], 
where = (1 + /i)/(l — /x). For L = 7.4, the symmetric (/i = 0, A = 1) minimum has lower 
energy than the asymmetric minima. For L = 8.4, the asymmetric minima are lower. 
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(c) $ (d) Energies for A; = 1, 2, 3 

Figure 4: Skyrmion properties for the system with V — (03)^[1 — (^s)^], as a function of L. 
In (a), the dot-dashed hne is the Bogomol'nyi bound; the sohd curve is the energy of the 
symmetric solution; and the dashed curve is the energy of the asymmetric solution. Figures 
(b) and (c) plot the quantities ^ and $ for the minimal-energy k — 1 solution. Figure 
(d) gives the energy for A; = 1 (solid curve), k — 2 (dashed curve) and k — 3 (dot-dashed 
curve), for the minimal-energy solutions. 
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